Current-induced spin polarization in spin-orbit-coupled electron systems 
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Current-induced spin polarization (CISP) is rederived in ballistic spin-orbit-coupled electron systems, based 
on equilibrium statistical mechanics. A simple and useful picture is correspondingly proposed to help under- 
stand the CISP and predict the polarization direction. Nonequilibrium Landauer-Keldysh formalism is applied to 
demonstrate the validity of the statistical picture, taking the linear Rashba-Dresselhaus [001] two-dimensional 
system as a specific example. Spin densities induced by the CISP in semiconductor heterostructures and in 
metallic surface states are compared, showing that the CISP increases with the spin splitting strength and hence 
suggesting that the CISP should be more observable on metal and semimetal surfaces due to the discovered 
strong Rashba splitting. An application of the CISP designed to generate a spin-Hall pattern in the inplane, 
instead of the out-of-plane, component is also proposed. 

PACS numbers: 72.25.Pn, 71.70.Ej, 85.75.-d 



I. INTRODUCTION 

The aim of preparing and controlling spins in all-electrical 
nonmagnetic devices has been shown to be possible in 
semiconducting bulk and two-dimensional electron systems 
(2DESs). ■ Besides the optical spin injection, a much more 
natural way of spin orientation is to make use of the spin-orbit 
(SO) coupling due to the lack of inversion symmetry of the un- 
derlying material. When passing an unpolarized electric cur- 
rent (electrons carrying random spins) through an SO-coupled 
material, spin-dependent consequences arise, among which 
two famous phenomena are the spin-Hall effect (SHE) 4 ' 5 ' 6 ' 7,8 ' 9 
and the current-induced spin polarization (CISP). 

In the CISP phenomenon, unpolarized electric current is ex- 
pected to be spin-polarized when flowing in a SO-coupled 
sample. This effect was first theoretically proposed in the 
early 90s. Edelstein employed linear-response theory to cal- 
culate the spin polarization due to an electric current in the 
presence of SO coupling linear in momentum, taking into 
account low-concentration impurities. Aronov and Lyanda- 
Geller 1 1 solved the quantum Liouville's theorem for the spin 
density matrix to show the CISP, taking into account scatter- 
ing as well. Recently, the CISP phenomenon has been exper- 
imentally proven. 12 ' 13 ' 14 Moreover, both the SHE and CISP 
have been observed at room temperature. 15 

In this paper we propose another viewpoint based on 
equilibrium statistical mechanics to explain the CISP in 
the absence of impurity scattering, for both bulk and two- 
dimensional systems. We show that the canonical ensemble 
average (CEA) of electrons moving with a wave vector k im- 
mediately prescribes a spin polarization antiparallel to the ef- 
fective magnetic field B e ff(k) stemming from the underlying 
SO coupling not necessarily linear in k, and hence explains the 
CISP. Correspondingly, a much simpler picture, compared to 
the early theoretical works of Refs. 10 and 1 1, helps provide 
a qualitative and straightforward explanation for the CISP: In 
an SO coupled 2DES without external magnetic field, an en- 
semble of rest electrons is unpolarized, while it becomes spin- 
polarized antiparallel to B e ff(k) when moving along k (see 
Fig. 1). 

To demonstrate the validness of this elementary statis- 



tical argument, spin and charge transports in finite-size 
four-terminal conducting 2DESs with Rashba and linear 
Dresselhaus [001] SO couplings, are numerically analyzed 
using the more sophisticated Landauer-Keldysh formalism 
(LKF), 16 ' 17 ' 18 allowing for nonequilibrium statistics. Good 
agreement between the analytical CEA and the numerical 
LKF will be seen, consolidating our statistical picture. In ad- 
dition to the semiconducting heterostructures, we also extend 
the analysis of the CISP to metal and semimetal surfaces, and 
compare the polarization strengths. Finally, an application of 
the CISP, resembling an inplane SHE, will be subsequently 
proposed. Throughout this paper, all the band parameters 
used in the LKF are extracted from experiments by match- 
ing the band structures calculated by the tight-binding model 
(and hence the density of states calculated by the LKF) with 
the experimentally measured ones. 19 

This paper is organized as follows. In Sec. II, we discuss 
the general properties of the system with SO coupling and de- 
rive the CISP in the ballistic limit using statistical mechanics. 
In Sec. Ill the LKF is applied partly to examine the validity 
of the statistical picture of the CISP introduced in Sec. II, and 
partly for further investigation. Summary of the present work 
will be given in Sec. IV. 
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FIG. 1: (Color online) Statistical picture of the current-induced spin 
polarization phenomenon. 
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II. ANALYTICAL DERIVATIONS 

Consider a SO-coupled system, subject to the single- 
particle Hamiltonian 

fi 2 k 2 

n = — i + s-n(k), (i) 

2m 

where m is the effective mass, 1 is the 2x2 identity ma- 
trix, S = (H/2)a is the spin operator, a = (cr x ,a y ,<r z ) 
being the Pauli matrix vector, and fi (k) = (e/mc)B e ff (k) 
is the momentum-dependent Larmor frequency vector, with 
B e ff (k) being the effective magnetic field stemming from the 
SO coupling. 20 



[001] model Hamiltonian. • With larger fc the second term 
in Eq. (4) — the fc 3 term — becomes important. We will come 
back to this later. For other quantum wells such as [110] and 
[111], the k vector given by Eq. (3) can be recast into a form 
that depends on the growth direction n of the 2DES. 24 (See 
also Ref. 20.) 

When writing the Larmor frequency vector as 

2a 

n R = —(kxn), (7) 

the linear Rashba model Hamiltonian 3 ' 20 ' 25 Hr = S • CIr = 
a(k x h) is recovered. Here a is the Rashba SO coupling 
constant. 



A. Larmor frequency vectors 



B. Time-reversal symmetry 



For III- V (zinc blende) bulk semiconductors," 1 the Larmor 
frequency in Eq. (1) is written as 22 



Q(k) 



rjh 2 



K. 



(2) 



where 77 is a dimensionless parameter specifying the spin-orbit 
coupling strength, E g is the band gap, and k is given by 



k x (k^ fc 2 ) 
ky(k z — k x ) 
kz(k x — k ) 



(3) 



Here fc^'s are the wave vector components along the crystal 
principle axes. 

When restricted to two-dimension, the component of the 
wave vector normal to the 2DES is averaged. For [001] quan- 
tum wells, one has fc 2 — > (fc 2 ) and k z — > (k z ) = (id z ) = to 
rewrite Eq. (3) as = [k x (k 2 y - (fc 2 )), fc y ((fc 2 ) - fc 2 ), 0], 

so that the Larmor frequency (2) takes the form 



^[001] _ 2 /3 
where (3 is defined by 



2d 

% (- kxi k y ,0) + (k X kl~kyklO) , (4) 



v h 2 



(2m?E g 



172<A£>=7<*£> 



(5) 



and is referred to as the Dresselhaus SO coupling constant. 
The 7 parameter (corresponding to bff c of Ref. 3) is 

material-dependent and is roug hly 27 eVA for both GaAs 
and InAs. 3 ' 

The first term in Eq. (4), 



fljp ' — — ( k x , k y , 0) . 



(6) 



is the linear Dresselhaus [001] term, which will dominate for 
small fc region. The corresponding SO term 7i^ 01 ' = S • 
fi^ 01 ' = (3(~k x a x +k y (T v ) is known as the linear Dresselhaus 



Before deriving the CISR we provide the following two in- 
trinsic properties of the Hamiltonian (1). First, we show that 
the contribution to the SO terms in solid is odd in fc due to 
time -reversal symmetry, which is also remarked in Ref. 3. 
For spin-1/2 systems subject to Hamiltonian (1), the energy 
dispersion can be written as 



E a (k) =E +aA k , 



(8) 



where Eq = h 2 k 2 /2m is the kinetic energy, a = ±1 is 
the spin state label, and Ak is the spin splitting due to SO 
coupling. In the absence of external magnetic field, the 
time-reversal symmetry is preserved, resulting in E + (k) = 
£L(-k),or, 



-A, 



(9) 



which implies that nonvanishing spin splitting A^ is odd in fc. 
Note that Eq. (9) also implies 



fi(-k) = -n(k). 



(10) 



which agrees with our intuition. Apparently, Eq. (10) is 
obeyed by all the previously reviewed Larmor frequency vec- 
tors. 

Second, we show (±, k|cr|±, k) = — (^f, k|<r|=F, k), where 
I £7, k) is the eigenstate of Hamiltonian (1). We begin with the 
Schrodinger equation, 



n\a,k) = ( 



h 2 k 2 

2m 



S-Q(k))|cr,k) = ^(k)|cr,k). (11) 



Comparing Eq. (1 1) with Eq. (8), we deduce S-0 (k) |cr, k) = 
crAk|cr, k), or, 

(tr,k|S-fi(k)|o-,k)=<rA k , (12) 
where \a, k) is assumed normalized. This implies 

(+,k|s-n(k)|+,k) = -(-, k|s-n(k)|- k). (13) 

Factoring out and canceling 51 (k) on both sides, we arrive at 
(+,k|<?|+,k> = -(- k|<?|-,k>. (14) 
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Equation (14) is a general property of Eq. (1) and is valid for 
systems with dispersions E a (k) = Eq + oAk, where the spin 
splitting A k is not necessarily linear in k. This property (14) 
will play a tricky role in the coming derivation of the CISP 
based on statistical mechanics in Sec. ??. 

Note that Eq. (14) is also a consequence of time-reversal 
symmetry (9), as one can easily prove as follows. Using Eq. 
(12) we rewrite Eq. (9) as 

(+,k|S-fi(k)|+,k) = (_ _k|S-fi(-k)|- -k). (15) 

Equation ( 1 2) also implies 

(cr,k|S -Q(k) |cr, k) = -{-er,k|S -0(k)| - cr,k) (16) 

when one regards oAk as — (— a) Ak. In addition, Eq. (9) 
implies 

(a, k|S ■ fi (k) |o-, k) = (cr, k S • n (k) |ff, -k) (17) 
because of 

(cr, k|S • (1 (k) |cr, k) = crA k = -crA_ k 

= (-ff, -k|S • Q, (-k) | -ff,-k) 

= -(ff,-k|S-0(-k)|cr,-k) (18a) 

= (ff,-k|S-0(k)|ff,-k), (18b) 

where Eqs. (16) and (10) are used in (18a) and (18b), respec- 
tively. Substituting Eqs. (10) and (17) into Eq. (15), we obtain 
Eq. (13), and hence the property (14). 



C. Current-induced spin polarization by canonical ensemble 
average 

Having seen the general properties of the Hamiltonian (1) 
under the time-reversal symmetry, we now derive the equilib- 
rium statistics version of the CISP. In quantum statistics, any 
physical quantity, say A, is expressed in terms of the quan- 
tum statistical average [A] — Tr(pA). Adopting the canonical 
ensemble, the average reads 



[A] 



Tr (e-«/ fe « T A) 



(19) 



where k B is the Boltzmann constant, T is temperature, v is 
a quantum number labeling the states, and E v is the eigenen- 
ergy of state v solved from Hamiltonian H. 

Now consider an unpolarized electron ensemble in a 2DES, 
subject to Hamiltonian (1). Our main interest here is the CEA 
of the spin operators of an ensemble of electrons, subject to an 
identical wave vector k. By this we mean that the summation 
in Eq. (19) runs over the spin index a only. This gives 



H Tr{e- n / ksT a) 
2J2a=± e- B -( k )/ fe B T ' 



Choosing the basis [cr, k) for the trace, one is led to 

7i£ CT e--M k y feBT (ff,k|ff|ff, k ) 



[S] k = 7T 



Using the property (14) and factoring out e fi2fc2 l 2mk BT f rom 
e -E a jk B T ^ we arr j ve at tne g enera i expression 

[S] k = -^tanh0;<+,k|ff|+,k>. (20) 

To re-express Eq. (20) in terms of the effective magnetic field 
B e ff (k), defined by 



B eff (k) EE — fi(k) = 2 

e 



fi(k) 



PB 



(21) 



we rewrite Eq. (12) with a = +1 as 



A, 



,k| ( ?|+,k)-B eff (k) = ^. 

Ms 



(22) 



Noting |(ff,k|ff|ff,k)| = 1 (unit vector) and |(7i/2)fi(k)| 
Ak, Eq. (22) implies 



<+,k|a|+,k> =B eff (k), 



(23) 



i.e., the direction of the effective magnetic field. Therefore, 
Eq. (20) can be written as 



[S] k = -^tanh-^B eff (k), 
2 k B T 



(24) 



^ CT e-B^(k)/fc B T 



which is exactly the analog of the CEA of electron spin in 
vacuum subject to an applied magnetic field. 26 

Equation (24) now has a transparent meaning: In the pres- 
ence of SO coupling, an ensemble of rest electrons (k — ► 0) is 
unpolarized since Ak^o = 0, while it becomes spin-polarized 
antiparallel to B e ff(k) when moving along k. This picture 
is schematically shown in Fig. 1. Moreover, the hyperbolic 
tangent factor tant^Ak/fesT 1 ) clearly predicts the decrease 
with T and the increase with Ak in the polarization magni- 
tude, and therefore explains two signatures of the CISP qual- 
itatively: (i) The CISP may persist up to the room temper- 
ature. Taking Ak ~ 3.68 meV from Ref. 14, one has 
tanh[A k /(fc B x 300K)]/tanh[A k /(fc s x 10 K)] « 14%. 
(ii) As (k) oc Vq (Ref. 13) implies Ak oc Vq, the magni- 
tude of the CISP governed by tanh(Ak/fcsT) is supposed to 
increase with the bias, as is experimentally proven. 12 



D. Explicit forms of current-induced spin polarization 

From Eq. (24), it is now clear that the direction of the CISP 
is given by the effective magnetic field direction B e ff(k). Al- 
ternatively, one can use the direction of the Larmor frequency 
vector, O(k), to describe the CISP direction since B e ff(k) and 
f2(k) are, by definition of Eq. (21), collinear. Therefore, the 
CISP direction in III-V bulk semiconductors is given by Eq. 
(3). 

For 2DES grown along [001] with Dresselhaus terms up to 
the fc 3 , Eq. (4) describes the effective magnetic field shown as 
Fig. 2, which simulates a 100-A-thick InGaAs quantum well 
with (k 2 z ) = 3.6 x 10~ 4 A~ 2 (Ref. 3). The CISP direction 
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k x (A" 1 ) 

FIG. 2: (Color online) Effective magnetic field of a 100-A-thick 
[001] InGaAs quantum well with {k 2 z ) = 3.6 x 10~ 4 ~ 2 . 




Rashba channel Linear Dresselhaus 

[001] channel 



FIG. 3: (Color online) Spin orientation in a 30a x 10a channel with 
a = 1 nm. Channels with linear Rashba model are considered in (a) 
and (c) while those with linear Dresselhaus [001] model are in (b) 
and (d). The direction of each sharp triangle represents the inplane 
spin vector (S)|| = ({S x ), (Sy)) of the local spin density. The size of 
the triangle depicts the magnitude of (S) m . Effective magnetic fields 
due to individually the Rashba and the Dresselhaus [001] fields are 
shown in the insets. 



is opposite to the effective magnetic field. Note that in Fig. 
2, the field distribution near the central region (small k) is 
dominated by the linear term (6) (cf. the right inset of Fig. 3). 

In the rest of this paper, we focus on the Rashba and linear 
Dresselhaus [001] terms. For effects with full SO terms in 
the Rashba-Dresselhaus systems, see Refs. 3 and 27. The 
composite Larmor frequency vector can be obtained by adding 
Eq. (6) with h = (0, 0, 1) and Eq. (7) together, 

n^ = a R (n = z) + nT ] 

= | [a (ky, -Ax, 0) + 13 (-k x , k v ,0)} . (25) 

The spin splitting linear in k takes the form At = \C\k 
with C = iae-W + pe**. Thus the CISP in linear Rashba- 
Dresselhaus [001] 2DESs is explicitly given by 

[S] k = ~2 k^f RD (26) 

E. Remark on effective mass 

In general, the inplane effective mass m of the electrons 
is not constant but depend strongly on k for realistic semi- 
conductor systems. However, in the long-wavelength limit 
kpa <c 1 (kp and a the Fermi wave vector and lattice con- 
stant, respectively), the effective mass, defined by the inverse 
of the second derivative of E(k)/h 2 with respect to k, is 
a constant due to the parabolic nature of E(k) solved from 
Hamiltonian (1). In this limit, even though the band structure 
can be anisotropic due to the interplay between different SO 
couplings (such as Rashba plus linear Dresselhaus [001]), the 
effective mass remains constant. In the present analysis, we 



work in this fcfa < 1 limit, within which the Hamiltonian (1) 
is valid. Interestingly, our CEA formulas such as Eq. (24) do 
not contain the dependence of m. 

Away from kpa <C 1 region, the energy dispersion is no 
longer parabolic, and the free-electron-like model Hamilto- 
nian (1) and hence the follow-up derivations fail. Analysis 
of the CISP phenomenon requires other formalisms such as 
the LKF, to be employed in the coming section. Nevertheless, 
we will not look further into the influence of the k-dependent 
effective mass on the CISP. 



III. NUMERICAL RESULTS: LANDAUER-KELDYSH 
FORMALISM 

To inspect the validity of the previously proposed statistical 
picture and further examine the CISP, we now perform local 
spin-density calculation in finite-size 2DESs attached to four 
normal metal leads by using the LKF. 16,17 ' 18 



A. Local spin densities in extreme Rashba and Dresselhaus 
[001] cases 

As a preliminary demonstration, Fig. 3 shows the position- 
dependent in-plane spin vectors (S)|j = ((S X ) T , (S y ) r ), with 
the local spin densities (S X ) T and (S y ) r calculated by the 
LKF. Here we adopt the finite difference method and dis- 
cretize the 30a x 10a channel, made of InGaAs/InAlAs 
heterostructure grown along [001], into a square lattice with 
lattice spacing a = 1 nm. Accordingly, this gives the kinetic 
and Rashba hopping strengths t = ti 2 /2ma 2 = 0.762 cV 
and t R = a /2a = 3.6 mcV, respectively. For the Dressel- 
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FIG. 4: (Color online) Local spin densities by LKF in a square 
Rashba-Dresselhaus [001] channel with (a) left to right, (b) left- 
bottom to right-top, (c) bottom to top, and (d) right-bottom to left-top 
bias configurations. Bias regime belongs to low: eVb = 2 meV. In- 
set: [S]^ vs. B^ 1 ' (k) in the k x -k y coordinate. 



haus SO coupling, we again assume the quantum well thick- 
ness d = 100 A and (k 2 z ) w (ir/d) 2 , and use 7 w 27cVA' 3 to 
give [see Eq. (5)] /3 = 7 (fc 2 ) « 2. 66 x 10~ 2 cV A, resulting 
in the Dresselhaus hopping strength to = /3/2a = 1.33 meV. 

Let us first consider the extreme cases, pure Rashba and 
pure Dresselhaus [001] channels. As expected, the spin vec- 
tors are mostly oriented antiparallel to B e ff(k), which is, for 
k || x, pointing to — y in the Rashba channel [Fig. 3(a)/(c) 
with low/high bias], and —x in the Dresselhaus [001] channel 
[Figs. 3 (b)/(d) with low/high bias]. Here (and hereafter) the 
low and high biases mean eVb = 2 meV and 0.2 cV, respec- 
tively, and we label the applied potential energy of ±eVo/2 as 
"±", and eVb = as "0" on each lead. Note that the spin dis- 
tribution, modulated by the charge distribution, forms stand- 
ing waves in the low bias regime since the electrons behaves 
quantum mechanically, while that in the high bias regime, i.e., 
the nonequilibrium transport regime, decays with distance. 
The polarization in the latter (high bias) is about two orders of 
magnitude stronger than the former (low bias). 

B. Consistency check: Analytical canonical ensemble average 
vs numerical Landauer-Keldysh formalism 

We now consider a four-lead square channel with coexist- 
ing Rashba and linear Dresselhaus [001] terms. The coupling 
constants are set identical to those introduced previously. Re- 
moving the four corner sites to avoid short circuit, the sam- 
ple size is (10 x 10 — 4)a 2 . To see if the CISP direction 
follows the opposite effective magnetic field for all k direc- 



tions, we change the current direction by applying different 
bias configurations. As shown in Figs. 4(a), (b), (c), and (d), 
the electrons flow from left to right, from left bottom to right 
top, from bottom to top, and from right bottom to left top, re- 
spectively. Other current directions are done in a similar way, 
but not explicitly shown here. In averaging the in-plane lo- 
cal spin densities (S x ) r and (S y ) r over all the lattice points 
at r within the conducting sample, we compare in the inset of 
Fig. 4 [S] k = ((S x ), (S y )) with the effective magnetic field 
B^ 1] (k) = (h/2fi B )nf D 1] (k), where ^ 1] (k) is given by 
Eq. (25). As expected by our statistical picture introduced 
in Sec. ??, [S]^ 1 ^ arrows are all opposite to O,^^ (k) for all 
k directions, despite some indistinguishably tiny differences. 
Note that the additive and destructive effects between the two 
SO terms are also observed at ±[110] and ±[110], respec- 
tively. Along ±[1 10] (±[110]), strongest (weakest) spin split- 
ting Ak, and hence the CISP magnitude [Eq. (24)], occur. 
Note that here we apply low bias. With high bias the results 
also agree perfectly with the CEA picture (not shown). 

C. Bias dependence of current-induced spin polarization 

Having shown that the statistical argument indeed works 
well, we next examine the bias dependence of the CISP, which 
is expected to be a proportional relation, as has been exper- 
imentally observed. 12 We return to Rashba channels. Spin 
densities, i.e., the total spin divided by the total area of 
the conducting channel, obtained via ^2 r (S y ) T / (Na 2 ) here 
with N being the number of total lattice sites in the con- 
ducting sample, are reported in Fig. 5 for sample widths 
W = 10a, 20a, 30a. Sample length is set L = 30a. Con- 
sistent to the experiment, the calculated spin densities in- 
crease with eVb- In addition, linear response within eVo < 
O.lto = 0.076 eV is clearly observed in all cases. Nonlin- 
earity enters when eVb grows so that nonequilibrium statistics 
dominates. Note that the calculated local spin density distri- 
bution satisfies the usual SHE symmetry, 18 so that we have 

£r<Sx) r = = and I E r (S)r| = Er(^>r. 

D. Comparison of current-induced spin polarization in 
semiconductor heterostructures and metal/semimetal surface 

states 

Next we extend the calculation of the spin density due to 
the CISP to other materials. In addition to semiconductor het- 
erostructures, 2DESs have been shown to exist also on metal 
surfaces supported by the surface states. Due to the loss of 
inversion symmetry, the metallic surfaces may exhibit Rashba 
spin splitting as well. 30,31 Here we consider three samples: 
54 x 18 nm 2 InGaAs/InAlAs heterostructure, 14.7 x 4.9 nm 2 
Au(l 1 1) surface, and 16.2 x 5.4 nm 2 Bi(l 1 1) surface. We ar- 
range the lead configuration of all the three samples as those 
in Fig. 3 and apply high bias. The sizes we choose here are 
to maintain roughly the same lattice site number N < 1000 
and keep the length-width ratio ~ 3. Note that realistic lat- 
tice structure are considered for the surface states [hexagonal 
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FIG. 5: (Color online) Bias dependence of spin densities induced by 
the CISP in Rashba 2DESs. 



( 



for Au(l 1 1) and honeycomb for Bi(l 1 1) bilayer], while finite- 
difference method based on the long-wavelength limit for the 
heterostructure is adopted. For introductory reviews of those 
surfaces, see Ref. 32 for noble metal surfaces, including gold, 
and Ref. for bismuth surfaces. 

Band parameters extracted from experiments and the spin 
densities calculated by the LKF are summarized in Table I. 
Clearly, the CISP increases with the Rashba parameter a. This 
suggests that the CISP (and actually also the SHE) should be 
more observable on these surfaces. The recently discovered 
Bi/Ag(lll) surface alloy that exhibits a giant spin splitting 
is even more promising, but we do not perform calculation for 
this interesting material here. 



E. Application of current-induced spin polarization: 
Generation of in-plane spin-Hall pattern 

Finally, we propose an experimental setup, as an applica- 
tion of the CISP, to generate an antisymmetric edge spin ac- 
cumulation in the inplane component, i.e., an inplane spin- 
Hall pattern. For simplicity, let us consider a Rashba 2DES 
with the parameters for the LKF calculation taken the same 



Material 


InGaAs/InAlAs 

(heterostructure) 


Au(lll) 

(surface state) 


Bi(lll) 

(surface state) 


m/mo 


0.050 


0.251 


0.340 


a(eVA) 


0.072 


0.356 


0.829 


E F - E b (eV) 


0.108 


0.417 


0.083 


Reference 


28 


30 


34 


CISP (l(T :i nrrr 2 ) 


0.240 


2.742 


8.382 



TABLE I: Summary of effective mass ratio m/mo, Rashba constant 
a, Fermi energy Ef (relative to the band bottom Eb), and the calcu- 
lated spin density due to CISP, for a set of materials. 



FIG. 6: (Color online) Mapping of the (a) local charge current den- 
sity, local spin densities (b) (S z ), (c) (S x ), and (d) (S y ) in a four- 
terminal square channel with a special bias arrangement. Unit in 
(b)-(c) is ti/2. 

as those in Fig. 5. Sample size is about 30 x 30 nm 2 . We 
apply high bias of eVo = 0.2 cV and arrange a special bias 
configuration. 

As shown in Fig. 6(a), unpolarized electron currents are 
injected from the left and right leads and are guided to the 
top and bottom ones. Under such design, the spin accumula- 
tion in (S z ) r exhibits merely a vague pattern [see Fig. 6(b)]. 
Contrarily, the pattern of (S x ) r shows not only antisymmetric 
edge accumulation in the channel but also magnitude much 
stronger than the out-of -plane component [see Fig. 6(c)]. This 
pattern is reasonably expected by the CISP due to the opposite 
charge flows along ±y at the top and bottom edges, and hence 
resembles an inplane SHE. 

In determining (S y ) r , Fig. 6(d) does not show a rotated 
pattern from (S x ) r due to the nonequilibrium transport. In 
the nonequilibrium transport regime, a distance apart from the 
source leads is required to induce the CISP, and therefore no 
significant (S y ) r is observed near the source (left and right) 
leads. This can be seen by comparing the local spin density 
distributions in the low-bias and high-bias regimes shown in 
Figs. 3(a) and 3(b), and Figs. 3(c) and 3(d), respectively. 



IV. SUMMARY 

In conclusion, we have rederived the CISP due to SO cou- 
pling in the absence of impurity scattering based on equilib- 
rium statistical mechanics. Correspondingly, a simple picture 
(Fig. 1) valid for both bulk structures and 2DESs is proposed 
to help qualitatively explain the CISP. Our explanation for the 
spin polarization of the moving electron ensemble in solid 
due to effective magnetic field is an exact analog to that of 
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the rest electron ensemble in vacuum due to external mag- 
netic field. 2 " The picture is further tested to work well even 
in the regime of nonequilibrium transport in finite-size sam- 
ples, by employing the numerical LKF. Extending the spin 
density calculation from the semiconductor heterostructure to 
metal and semimetal surface states, our calculation confirms 
that the polarization increases with the SO coupling strength, 
and hence suggests that the CISP should be more observable 
on metal and semimetal surfaces with stronger Rashba SO 
coupling. 30, ,35 As an application of the CISP, we also suggest 
an interesting bias configuration for the four-terminal setup to 



generate inplane SHE [Fig. 6(c)]. 
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